Mapping structural variants to rare disease genes using long-read whole genome sequencing and trait-relevant polygenic scores.

Recent studies have revealed the pervasive landscape of rare structural variants (rSVs) present in human genomes. rSVs can have extreme effects on the expression of proximal genes and, in a rare disease context, have been implicated in patient cases where no diagnostic single nucleotide variant (SNV) was found. Approaches for integrating rSVs to date have focused on targeted approaches in known Mendelian rare disease genes. This approach is intractable for rare diseases with many causal loci or patients with complex, multi-phenotype syndromes. We hypothesized that integrating trait-relevant polygenic scores (PGS) would provide a substantial reduction in the number of candidate disease genes in which to assess rSV effects. We further implemented a method for ranking PGS genes to define a set of core/key genes where a rSV has the potential to exert relatively larger effects on disease risk. Among a subset of patients enrolled in the Genomic Answers for Kids (GA4K) rare disease program (N=497), we used PacBio HiFi long-read whole genome sequencing (lrWGS) to identify rSVs intersecting genes in trait-relevant PGSs. Illustrating our approach in Autism (N=54 cases), we identified 1,827 deletions, 158 duplications, 619 insertions, and 14 inversions overlapping putative core/key PGS genes. Additionally, by integrating genomic constraint annotations from gnomAD, we observed that rare duplications overlapping putative core/key PGS genes were frequently in higher constraint regions compared to controls (P = 2×10-04). This difference was not observed in the lowest-ranked gene set (P = 0.18). Overall, our study provides a framework for the annotation of long-read rSVs from lrWGS data and prioritization of disease-linked genomic regions for downstream functional validation of rSV impacts. To enable reuse by other researchers, we have made SV allele frequencies and gene associations freely available.

Genomic prediction of relapse in recipients of allogeneic haematopoietic stem cell transplantation.

Allogeneic haematopoietic stem cell transplantation currently represents the primary potentially curative treatment for cancers of the blood and bone marrow. While relapse occurs in approximately 30% of patients, few risk-modifying genetic variants have been identified. The present study evaluates the predictive potential of patient genetics on relapse risk in a genome-wide manner. We studied 151 graft recipients with HLA-matched sibling donors by sequencing the whole-exome, active immunoregulatory regions, and the full MHC region. To assess the predictive capability and contributions of SNPs and INDELs, we employed machine learning and a feature selection approach in a cross-validation framework to discover the most informative variants while controlling against overfitting. Our results show that germline genetic polymorphisms in patients entail a significant contribution to relapse risk, as judged by the predictive performance of the model (AUC = 0.72 [95% CI: 0.63-0.81]). Furthermore, the top contributing variants were predictive in two independent replication cohorts (n = 258 and n = 125) from the same population. The results can help elucidate relapse mechanisms and suggest novel therapeutic targets. A computational genomic model could provide a step toward individualized prognostic risk assessment, particularly when accompanied by other data modalities.

Bone mineral density, osteoporosis, and osteoporotic fractures: a genome-wide association study.

BACKGROUND: Osteoporosis is diagnosed by the measurement of bone mineral density, which is a highly heritable and multifactorial trait. We aimed to identify genetic loci that are associated with bone mineral density. METHODS: In this genome-wide association study, we identified the most promising of 314 075 single nucleotide polymorphisms (SNPs) in 2094 women in a UK study. We then tested these SNPs for replication in 6463 people from three other cohorts in western Europe. We also investigated allelic expression in lymphoblast cell lines. We tested the association between the replicated SNPs and osteoporotic fractures with data from two studies. FINDINGS: We identified genome-wide evidence for an association between bone mineral density and two SNPs (p<5x10(-8)). The SNPs were rs4355801, on chromosome 8, near to the TNFRSF11B (osteoprotegerin) gene, and rs3736228, on chromosome 11 in the LRP5 (lipoprotein-receptor-related protein) gene. A non-synonymous SNP in the LRP5 gene was associated with decreased bone mineral density (rs3736228, p=6.3x10(-12) for lumbar spine and p=1.9x10(-4) for femoral neck) and an increased risk of both osteoporotic fractures (odds ratio [OR] 1.3, 95% CI 1.09-1.52, p=0.002) and osteoporosis (OR 1.3, 1.08-1.63, p=0.008). Three SNPs near the TNFRSF11B gene were associated with decreased bone mineral density (top SNP, rs4355801: p=7.6x10(-10) for lumbar spine and p=3.3x10(-8) for femoral neck) and increased risk of osteoporosis (OR 1.2, 95% CI 1.01-1.42, p=0.038). For carriers of the risk allele at rs4355801, expression of TNFRSF11B in lymphoblast cell lines was halved (p=3.0x10(-6)). 1883 (22%) of 8557 people were at least heterozygous for these risk alleles, and these alleles had a cumulative association with bone mineral density (trend p=2.3x10(-17)). The presence of both risk alleles increased the risk of osteoporotic fractures (OR 1.3, 1.08-1.63, p=0.006) and this effect was independent of bone mineral density. INTERPRETATION: Two gene variants of key biological proteins increase the risk of osteoporosis and osteoporotic fracture. The combined effect of these risk alleles on fractures is similar to that of most well-replicated environmental risk factors, and they are present in more than one in five white people, suggesting a potential role in screening.

Association study between the CX3CR1 gene and asthma.

CX3CR1, a fractalkine receptor, mediates cell-adhesive and migratory functions in inflammation. Based on CX3CR1 expression observed in bronchial tissues of asthmatic subjects, we hypothesized that genetic variation at this locus may affect susceptibility to asthma. We carried out an association study and a haplotypic analysis with selected polymorphisms of the CX3CR1 in a familial asthmatic sample from a founder population. Genetic analyses performed by FBAT software showed five CX3CR1 single nucleotide polymorphisms (rs938203, rs2669849, rs1050592, T280M and V249I) with significant associations between their common alleles and asthma (P<0.004) in a dominant model. A haplotype formed with common alleles of rs1050592, T280M and V249I is also overtransmitted in asthmatic subjects (P=0.005) under a dominant model. The associations of V249I and rs2669849 have been validated in an independent case-control sample. For V249I, odds ratios (OR) are 2.16 (common homozygous) and 2.11 (heterozygous) in dominant model (P=0.031). For rs2669849, OR are 2.75 (common homozygous) and 1.86 (heterozygous) in additive model (P=0.007) and dominant model (P=0.059). These results suggest an asthma protective effect of the minor alleles in healthy control carriers. Further functional studies of CX3CR1 are needed to document its role in the pathophysiology of asthma.

Folate intake, plasma folate and homocysteine status in a random Finnish population.

OBJECTIVE: To assess the folate status of Finnish adults using plasma folate and homocysteine as biomarkers and to evaluate dietary and supplementary folate intakes. MATERIALS AND METHODS: Plasma folate, vitamin B(12) and total homocysteine (tHcy) were determined in a random sample of 643 subjects aged 25-74 y living in the Helsinki area. The methylenetetrahydrofolate reductase (MTHFR)-genotypes were analyzed from a subsample (n=394). Dietary intake data by 24 h recall and use of vitamin supplements were collected. RESULTS: Plasma folate was normal (>/=5 nmol/l) in 99% of subjects and optimal (>/=8 nmol/l) in terms of a minimum tHcy in 90%. Mean plasma folate of non-supplement users was 13.7 and 12.9 nmol/l and tHcy 11.3 and 9.2 micro mol/l for men and women, respectively. Elevated tHcy (>14 micro mol/l) was found in 11% of subjects. Homozygote frequency for MTHFR genotype TT was 5.0% and their plasma tHcy was 14.8 micro mol/l compared to the mean of the other subjects, 10.5 micro mol/l, P<0.05. The mean dietary folate intake was 241 micro g/day (29 micro g/MJ of energy) for men and 205 micro g/day (33 micro g/MJ) for women, respectively. The main dietary sources of folate were vegetables 12%, wholemeal ryebread 11%, fruits 10%, and potato 10%. Regular supplement users (n=97) received on average 207 micro g folic acid per day from supplements. CONCLUSIONS: The folate status of Finnish adults seems to be adequate according to energy adjusted folate intake, plasma folate and homocysteine. The MTHFR homozygote frequency was low compared to other countries. Regular use of supplementary folic acid less than 300 micro g increased plasma folate, but supplemental folic acid over 300 micro g was required to lower tHcy values significantly.

Chromosome 19q13 and multiple sclerosis susceptibility in Finland: a linkage and two-stage association study.

Several studies have previously provided some albeit weak evidence for linkage or association between chromosome 19q13 and multiple sclerosis (MS) susceptibility. We performed a two-stage association analysis with 19 markers spanning 7 Mb/5.5 cM of 19q13. In stage 1 analysis (135 MS families) allelic and haplotypic associations were found with markers within or close to the ApoE-ApoC subregion. These observations were taken as a hypothesis, which was tested in stage 2 in 125 families. However, none of the initial associations were replicated suggesting that they were most likely due to chance. Linkage analysis was performed in 27 Finnish multiplex families using 10 microsatellites spanning 23 Mb/24 cM of 19q13. DNA was available from 72 MS patients and 150 unaffected relatives. Parametric and non-parametric linkage analyses did not provide evidence for linkage when all families were tested. After stratifying the families according to HLA-DR15 there was weak evidence for linkage to the 19q13.1 subregion in DR15 negative families (LOD(max)=1.8). Taken together these results do not support a major role of chromosome 19q13.2-q13.3 in MS susceptibility among Finnish MS patients, whereas conclusions on the 19q13.1 subregion are less clear and this region requires further study.

Dissecting a population genome for targeted screening of disease mutations.

Compared to mixed populations, population isolates such as Finland show distinct differences in the prevalence of disease mutations. However, little information exists of the differences on the prevalence of different disease alleles in regional populations with different history of multiple bottlenecks. We constructed a DNA-array and monitored the prevalence of 31 rare and common disease mutations underlying 27 clinical phenotypes in a large population-based study sample. Over 64 000 genotypes were assigned in 2151 samples from four geographical areas representing early and late settlement regions of Finland. Each sample was analyzed in duplicate and a total of 142 000 array-derived genotyping calls were made. On average one in three individuals was found to be a carrier of one of the 31 monitored mutations. This should remove fears of the stigmatizing effect of a carrier-screening program monitoring multiple diseases. Regional differences were found in the prevalence of mutations, providing molecular evidence for the deviating population histories of regional subisolates. The mutations introduced early into the population revealed relatively even distribution in different subregions. More recently introduced rare mutations showed local clustering of disease alleles, indicating the persistence of population subisolates and the effect of multiple bottlenecks in molding the population gene pool. Regional differences were observed also for common disease alleles. Such precise information of the carrier frequencies could form the basis for targeted genetic screens in this population. Our approach describes a general paradigm for large-scale carrier-screening programs also in other populations.

Y-chromosomal SNPs in Finno-Ugric-speaking populations analyzed by minisequencing on microarrays.

An increasing number of single nucleotide polymorphisms (SNPs) on the Y chromosome are being identified. To utilize the full potential of the SNP markers in population genetic studies, new genotyping methods with high throughput are required. We describe a microarray system based on the minisequencing single nucleotide primer extension principle for multiplex genotyping of Y-chromosomal SNP markers. The system was applied for screening a panel of 25 Y-chromosomal SNPs in a unique collection of samples representing five Finno--Ugric populations. The specific minisequencing reaction provides 5-fold to infinite discrimination between the Y-chromosomal genotypes, and the microarray format of the system allows parallel and simultaneous analysis of large numbers of SNPs and samples. In addition to the SNP markers, five Y-chromosomal microsatellite loci were typed. Altogether 10,000 genotypes were generated to assess the genetic diversity in these population samples. Six of the 25 SNP markers (M9, Tat, SRY10831, M17, M12, 92R7) were polymorphic in the analyzed populations, yielding six distinct SNP haplotypes. The microsatellite data were used to study the genetic structure of two major SNP haplotypes in the Finns and the Saami in more detail. We found that the most common haplotypes are shared between the Finns and the Saami, and that the SNP haplotypes show regional differences within the Finns and the Saami, which supports the hypothesis of two separate settlement waves to Finland.

A system for specific, high-throughput genotyping by allele-specific primer extension on microarrays.

This study describes a practical system that allows high-throughput genotyping of single nucleotide polymorphisms (SNPs) and detection of mutations by allele-specific extension on primer arrays. The method relies on the sequence-specific extension of two immobilized allele-specific primers that differ at their 3'-nucleotide defining the alleles, by a reverse transcriptase (RT) enzyme at optimized reaction conditions. We show the potential of this simple one-step procedure performed on spotted primer arrays of low redundancy by generating over 8000 genotypes for 40 mutations or SNPs. The genotypes formed three easily identifiable clusters and all known genotypes were assigned correctly. Higher degrees of multiplexing will be possible with this system as the power of discrimination between genotypes remained unaltered in the presence of over 100 amplicons in a single reaction. The enzyme-assisted reaction provides highly specific allele distinction, evidenced by its ability to detect minority sequence variants present in 5% of a sample at multiple sites. The assay format based on miniaturized reaction chambers at standard 384-well spacing on microscope slides carrying arrays with two primers per SNP for 80 samples results in low consumption of reagents and makes parallel analysis of a large number of samples convenient. In the assay one or two fluorescent nucleotide analogs are used as labels, and thus the genotyping results can be interpreted with presently available array scanners and software. The general accessibility, simple set-up, and the robust procedure of the array-based genotyping system described here will offer an easy way to increase the throughput of SNP typing in any molecular biology laboratory.

Array-based multiplex analysis of candidate genes reveals two independent and additive genetic risk factors for myocardial infarction in the Finnish population.

We analysed common variants of eight genes implicated previously as risk factors for coronary heart disease or myocardial infarction (MI) in a cross-sectional study on patients with a history of MI and in carefully matched controls from the Finnish population. The most common low density lipoprotein receptor mutations in Finland were also included in our analysis. Multiplex genotyping of the target genes was performed using a specific and efficient array-based minisequencing system. The 4G allele of the plasminogen activator inhibitor gene (P < 0.05) and the PlA2 allele of the glycoprotein IIIa gene (P < 0.01) were associated with an increased risk of MI in our study population. We analysed the combined effect of these risk alleles and found that the concurrent carrier status of the two genetic variants conferred a high risk for the development of MI in our sample (OR = 4.5, P = 0.001), which was particularly prominent in male subjects (OR = 6.4, P = 0.0005). This study demonstrates the application of a new powerful tool for genome analysis to yield information on the inherited determinants of susceptibility to MI. The observation of two separate genes contributing an additive risk of developing MI exemplifies the advantages of multiplex analysis of genetic variation.

Contribution of the CCR5 and MBL genes to susceptibility to HIV type 1 infection in the Finnish population.

Homozygosity for a variant of a chemokine receptor gene (CCR5) has been shown to protect from HIV-1 infection. Variants of the mannose-binding lectin (MBL) gene have been suggested to predispose to HIV-1 infection. These genetic variants and their possible role in susceptibility to HIV-1 infection were studied in sample material from more than 300 Finnish HIV-infected and control individuals. The genotyping was carried out efficiently using a novel, primer extension assay in a miniaturized oligonucleotide array format. Homozygotes for the protective deletion allele of the CCR5 gene were found only in the control group, and the frequency of the allele was high in the Finnish population. Homozygosity for the MBL variant alleles was enriched significantly in the HIV-1-infected group, thus providing further evidence for the harmfulness of MBL variant homozygosity in HIV-1 transmission.

Golli-MBP gene in multiple sclerosis susceptibility.

Multiple sclerosis (MS) is an oligo- or polygenic disease but no specific susceptibility genes have been identified so far. In the Finnish population we have previously found evidence for linkage between MS and the myelin basic protein gene (here called Golli-MBP gene) suggesting that either Golli-MBP or another gene in its vicinity contributes to MS suceptibility. Here we have screened the Golli-MBP gene for nucleotide variations and carried out multipoint association analyses in a Finnish case-control data-set as well as in an independent data-set composed of 151 MS families from Finland and Sweden. In both data-sets we found association between MS and alleles in the 1.27 kilobase (kb) range at a tetranucleotide repeat element (TGGA)n which is located 1 kb upstream of the MBP exon 1. Haplotype analyses suggested that the MS-associated 1.27 kb alleles can be split into predisposing and non-predisposing variants and provided evidence that the candidate DNA region contributing to MS susceptibility should be located at the Golli-MBP gene within a 20-25 kb segment that was conserved in the predisposing haplotypes. These findings suggest a role for the Golli-MBP locus in MS susceptibility, at least in a subset of patients, and serve as a basis for highly focused attempts to identify predisposing mutation(s).

Immune system genes in multiple sclerosis: genetic association and linkage analyses on TCR beta, IGH, IFN-gamma and IL-1ra/IL-1 beta loci.

The role of genetic factors in the etiology of multiple sclerosis (MS) has been clearly demonstrated but the loci determining susceptibility to this disease remain largely unidentified. A contribution from several immune system genes has been suggested based on animal models and association/linkage analyses on MS patients and families. With the exception of the findings from the HLA complex, studies on candidate immune system genes have provided controversial results. Here we have performed genetic association and linkage analyses on four chromosomal regions containing immune system genes. A possible role for each of these loci in MS has been previously suggested. In data-sets derived from the Finnish population we found no evidence for contribution of the T-cell receptor beta chain (TCR beta chromosome 7q35), immunoglobulin heavy chain (IGH chromosome 14q32), interferon-gamma (IFN-gamma chromosome 12q14-q15) or interleukin-1 receptor antagonist/interleukin-1 beta (IL-1ra/IL-1 beta chromosome 2q14-q21) loci in the genetic susceptibility to MS.

Minisequencing: a specific tool for DNA analysis and diagnostics on oligonucleotide arrays.

We describe a method for multiplex detection of mutations in which the solid-phase minisequencing principle is applied to an oligonucleotide array format. The mutations are detected by extending immobilized primers that anneal to their template sequences immediately adjacent to the mutant nucleotide positions with single labeled dideoxynucleoside triphosphates using a DNA polymerase. The arrays were prepared by coupling one primer per mutation to be detected on a small glass area. Genomic fragments spanning nine disease mutations, which were selected as targets for the assay, were amplified in multiplex PCR reactions and used as templates for the minisequencing reactions on the primer array. The genotypes of homozygous and heterozygous genomic DNA samples were unequivocally defined at each analyzed nucleotide position by the highly specific primer extension reaction. In a comparison to hybridization with immobilized allele-specific probes in the same assay format, the power of discrimination between homozygous and heterozygous genotypes was one order of magnitude higher using the minisequencing method. Therefore, single-nucleotide primer extension is a promising principle for future high-throughput mutation detection and genotyping using high density DNA-chip technology.

Multiplex, fluorescent, solid-phase minisequencing for efficient screening of DNA sequence variation.

We developed a multiplex, solid-phase minisequencing method to detect multiple single-nucleotide polymorphisms in an undivided sample. The amplified DNA templates are first captured on a manifold. Then, with multiple minisequencing primers of various sizes, single-nucleotide extension reactions are carried out simultaneously with fluorescently labeled dideoxynucleotides. The size of the extended product, determined by using a DNA sequencing instrument, defines the site of the polymorphisms, and the incorporated nucleotide gives the identity of the nucleotide at each site. HLA-DQA1 typing was used as a model system to evaluate the method. The DR2 subgroup of the HLA-DRB1 gene was typed along with the DQA1 gene to demonstrate the feasibility of the method in analyzing multiple genes at multiple sites simultaneously. The method is generally applicable for screening any single-nucleotide polymorphisms or point mutations, and its manifold format allows practical handling of large numbers of samples.